Electron-electron correlations in a dynamical impurity system with a Fermi edge 

singularity 
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We study spatial correlations in the ground state of a one-dimensional electron gas coupled to 
a dynamic quantum impurity. The system displays a non-trivial many-body effect known as the 
Fermi edge singularity: transitions between discrete internal states of the impurity have a power- 
law dependence on the internal energies of the impurity states. We present compact formulas for 
the static current-current correlator and the pair correlation function. These reveal that spatial 
correlations induced by the impurity decay slowly (as the third inverse power of distance) and have 
a power-law energy dependence, characteristic of the Fermi edge singularity. 
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I. INTRODUCTION 

The situation where a localized scatterer, with discreet 
internal quantum mechanical degrees of freedom, is cou- 
pled to an electron gas, is a common occurrence both in 
bulk condensed matter systemsi and in nano-meter scale 
electronic devices.—"— Many non-trivial features of these 
systems can be explained with the aid of simple dynamic 
quantum impurity models. '"^ For several such models, non- 
perturbative results have been obtained and these results 
have played a major role in our understanding of strongly 
correlated electron systems. 

In this context, a phenomenon known as the Fermi 
edge singularity^ has been very influential. Its essential 
ingredients are a Fermi gas that is initially in a stationary 
state and a local scattering potential that is abruptly 
switched on at a time tg. For times larger than to, the 
Fermi gas is no longer in a stationary state. To see just 
how far from stationary a typical initial state is, it is 
useful to consider the overlap between the ground states 
of the Hamiltonian before and after tQ. It is found that 
the overlap is zero, and this is known as the orthogonality 
catastrophcfi Since the initial state is far from stationary, 
the Fermi sea is severely shaken up;^ the local scattering 
potential creates a multitude of particle hole excitations. 

The original theory^iiS was formulated to account for 
singularities in the photo-emission and -absorption spec- 
tra in some metals. For this problem the natural quan- 
tities to investigate are impurity transition rates. These 
where found to have a power-law dependence on the in- 
ternal energies of the impurity levels. It was subsequently 
shown that the same physics applies to transport through 
a barrier containing an impurityii level, or through 
a quantum dotr^^^ In recent work, non-equilibrium 
set-ups in which the Fermi sphere is replaced by a 
non-equilibrium distribution has been investigated^^"— . 
Complications to the basic theory, for instance what hap- 
pens if the impurity induces resonant scattering, was also 
considered. Apart from impurity transition rates, cur- 
rent and noise^^ as well as the quench dynamics of the 
electron gas22^ have been investigated. 



These studies clearly indicate that the dynamical im- 
purity induces significant correlations in the electron gas. 
A system that displays the Fermi edge singularity can 
therefore be expected to possess a non-trivially corre- 
lated ground state. However, as far as we know, this has 
not yet been studied. Perhaps the reason is that some 
common realizations of the Fermi edge singularity rely on 
explicitly time dependent Hamiltonians. In this work we 
consider a realization with a time independent Hamilto- 
nian. This allows us to investigate ground state correla- 
tions. Correlations in time, such as current noisCfS? have 
already been studied. However, spatial correlations, a 
topic that is receiving significant attention in other im- 
purity systems)^"— have not yet been addressed. We 
therefore do so in this work. 

We consider a set-up in which a localized impurity in- 
teracts with one-dimensional electrons in a single chan- 
nel with a linear dispersion relationi^ An example of the 
type of impurity we have in mind, is the charge qubit 
formed when a single electron is trapped in a double 
quantum doti^i^ Our model can also be realized with a 
Joscphson charge qubit. "^^ Our main result is a compact 
formula for the current-current correlator {Sj{x)Sj{y)) 
(Eg. 14. 26]) . that is valid at large separations. The impu- 
rity induces no additional correlations between electrons 
on the same side of the impurity. For electrons on oppo- 
site sides of the impurity, impurity induced correlations 
decay as l/\x — y\^. Owing to the Fermi-edge singularity, 
the correlator has a power-law dependence on the inter- 
nal energies of the impurity. The correlations persist be- 
yond the range of the interaction between the electron gas 
and the impurity. This is in contrast to the expectation 
values of single particle observables such as the density 
{p{x)), which is unaffected by the impurity for x outside 
the range of the interaction. Because of the linear disper- 
sion relation, the current-current correlator is essentially 
equal to the pair correlation function {p{x)p{y)). This 
insight allows for a simple interpretation of our results. 
If an electron is detected at x to the left of the impurity, 
the likelihood of finding another at y to the right of the 
impurity is larger than in the absence of the impurity by 
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an amount ^ l/\x — y\^ and has a power-law dependence 
on the internal energies of the impurity, characteristic of 
the Fermi edge singularity. 

We present the work as follows: We define the model 
by means of its Hamiltonian (Sec. Ill Al ) Since we are 
dealing with one dimensional electrons with a linear dis- 
persion relation, it is convenient to bosonize them. In 
Sec. Ill BI we collect the relevant bosonization results. We 
are interested in the model because it displays a Fermi 
edge singularity. We start Sec. Ill CI by reviewing how 
the conventional theory of the Fermi edge singularity ap- 
plies to the transition rates between internal states of the 
impurity. We also take the first step towards studying 
ground state properties by considering the probability to 
find the impurity in a given internal state, if the sys- 
tem as a whole is in the ground state. These occupation 
probabilities turn out to play a central role in determin- 
ing spatial correlations between electrons in the presence 
of the impurity. We show that the Fermi edge singular- 
ity is also manifested in the occupation probabilities. As 
already stated, our main aim is to investigate spatial cor- 
relations among the electrons. But before doing so, we 
need to consider the expectation values of single parti- 
cle observables. In Sec. lIIII we therefore investigate the 
average density {p{x)). We show that the total number 
of particles displaced by the impurity obeys a general- 
ization of Friedel's sum rule. (A generalized Friedel sum 
rule has previously been shown to hold for another dy- 
namic impurity system, namely the Anderson impurity 
model.) Our results for the current-current correlator is 
presented in Sec. lIVI and for the pair correlation function 
in Sec.|Vl Sec. lVIl contains a self-contained summary of 
our main result and conclusions. 



II. MODEL 

Our model describes a localized impurity with an in- 
ternal state space spanned by the orthonormal basis 
B = {\a) |q! = 1, . . . , M}. It interacts with electrons in 
a single channel with a linear dispersion relation — k. 
The model is inspired by experimental setups in which 
elecrons in a coherent conductor are coupled to a charge 
qubit [M — 2).^^'^^ Since it does not substantially com- 
plicate our analysis, we consider arbitrary M . 

A. Hamiltonian 

The Hamiltonian for our model is H = Hq+Ht, where 
i/o-JI-^"® 1")^"!' HT = Y,7afi\a){P\, (2.1) 

a ap 
fL/2 

Ha= dx'ip\x)[-idx + Va{x)]ipix)+ea. (2.2) 

J-L/2 

Here Ha describes the fermions when the impurity is in 
state I a) and Sa is the internal energy of the impurity 



state la). The fermion creation and annihilation opera- 
tors ip'^ (x) and 'ipix) obey periodic boundary conditions 
on the interval (— i/2,i/2]. We will eventually send 
the system size L to infinity. The electrostatic potential 
Va{x) depends on the state of the impurity. Since our 
aim in this work is to investigate many-body correlations 
that persist beyond the range of the interaction between 
the impurity and the Fermi gas, it is convenient to as- 
sume that all the potentials Va have bounded support. 
We define the scattering region as the shortest interval 
I = [x-,x+] such that Va{x) = for all a and all x ^ I. 
The term Ht induces tunnelling between impurity states. 
Without loss of generality, we set jaa = 0. 

The ground state \Fa) of the fermion Hamiltonian Ha 
is a Fermi-sea in which all (the infinitely many) negative 
energy orbitals associated with the single particle Hamil- 
tonian —idx + Va{x) are occupied, and all positive energy 
orbitals are empty. 

The scattering states of —idx + Va{x) are of the form 
exp[ikx~i dx'va{x')]. This leads to an identification 
of 

1 f°° 

5a = -7: / dxva{x) (2.3) 
^ J -00 

as the scattering phase shift associated with the potential 
Va{x). Phase shifts play an important role in impurity 
physics. 

B. Bosonization 

The infinite Fermi sea represented by a state such 
as \Fa) must be handled with care. Finite observable 
quantities are often represented as the differnece be- 
tween infinite quantities. In such calculations, an un- 
controlled rearrangement of terms can lead to incorrect 
answers .^^1'^^ For one-dimensional Fermi systems, such 
as the one we are studying, a standard and elegant 
method to treat the infinitely deep Fermi sea correctly, 
is provided by bosonization.>24i24 In fact, one of the pi- 
oneering applications of bosonization was to the Fermi 
edge singularityj^ Through bosonization, our model is 
mapped onto a spin-boson model^ with a bath that is 
Ohmic at low energies For our purposes the following 
bosonization results are required. The Fourier compo- 
nents of the density operator 
rL/2 

Pg^ / dxe-""=i;'<{x)i^{x), (2.4) 

J-L/2 

with q = 21171/ L, n — 0, ±1, ±2, . . . obey bosonic com- 
mutation relations 

[Pq^ Pq'] ^ TT^q.-q' ■ (2-5) 

zvr 

The fermion Hamiltonians Ha can be expressed in terms 
of pq as Ha — ha + Eoa whcrc 

'.-i:|(^-.+^)(^.+l')^ f^.^) 

q>0 
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Here Va,q = J^j^^^^'^^^ ^'^^Va{x) are the Fourier compo- 
nents of the potentials Va {x) and 



Oq 



L 



E 

q>0 



2tt 



2tt 



(2.7) 



is the ground state energy of Ha as a function of po, 
the number of electrons in the system, up to an a- 
independent constant. We will work with a fixed number 
of particles. The number operator po, and therfore also 
, can then be treated as real numbers in stead of op- 
erators. Note that Eqc can be varied independently of 
ha by adjusting the internal energy of the impurity 
state a. We define the energy differences 



For g > 0, 



i^af! — Eoa — Eop. 



P,\F.)^-^\Fa) 



(2.8) 



(2.9) 



i.e. the Fermi sea \Fa) is an eigcnstate (coherent state) 
of the bosonic annihilation operator pq. 



C. Fermi edge singularity 

The Fermi edge singularity manifests as power law sin- 
gularities in inelastic transition rates between internal 
states of the impurity, as a function the energy differences 
u>ap- Denote the many-body eigenstates of 77^ by 
and let Em/s be the associated energies. Assuming an ini- 
tial state \Fa) (8) \a), the total transition rate W^a(a;^a) 
to configurations with the impurity in state calcu- 
lated to second order in Ht using Fermi's golden rule, is 
given by 

Wpail^Pa) 

= 2n\jfia\^Y.S{E^p - Eoa)\ {Fa\ mp) p 

m 

/oo 
dte^(^-^~^o.)t^{Fa\m,)f 
-OO 



{Fa 



^ihf,t^-iha,t 



\Fa). (2.10) 



The canonical theory of the Fermi edge singularity pro- 
vides an expression (the so-called closed loop factor) 
for the expectation value (i^^l e^'^'^'e"*''"* The 
closed loop factor is analytic in the upper half of the 
complex plain. For large \t\, its asymptotic behaviour is 

of the form (iM) "i^ where 

Aap = -{Sa-5p)/n. (2.11) 

A branch cut in the lower half of the complex t plane is 
understood, and the branch with arg(t) = for real pos- 
itive t is implied. A is a large energy scale. In our model. 



with an infinitely deep Fermi-sea, it is of the order of the 
inverse of the range of the interaction of the impurity and 
the Fermi gas. The interpretation is that an interaction 
with range A~^ can excite particle-hole pairs with ener- 
gies ~ A. In models where the inverse of the interaction 
range is larger than the Fermi energy, measured from the 
bottom of the Fermi sea, A is of the order of the Fermi 
energy. From the asymptotics and analiticity structure 
of the closed loop factor, the asymptotic behaviour of the 
rates Wpa can be extracted at small energy differences 



Wap{iOl3a) = |7q;3|^C'q^A 



1 



(2.12) 

Cap is a constant that tends to a value of order one when 
Aq.^ — 0. The result for the rate breaks down at energies 
ujap corresponding to excitations with wave-lengths w"^ 
short enough to resolve the spatial structure of the in- 
teraction potentials Va(x) and vp(x). For energies larger 
than A, the transition rate becomes exponentially sup- 
pressed. There is also another constraint. For perturba- 
tion theory in ^ap to be valid, the ratio Wi3a(!-0(3a)/<^ai3 
must be much smaller than onei^ For A^^ < 2 this im- 
poses the constrainli^ that 



> 



Ca 



A2 



(2.13) 



The point ujpa = is known as the Fermi edge threshold. 
More information regarding the above statements can be 
found in Ref.[2^, where we analysed transition rates for 
essentially the same model in detail. The fact that, for 
A^Q < 1, the transition rate W^^^^ becomes large close to 
the threshold, while for A^q, > 1 it becomes small close 
to the threshold, is an instance of what is known as a 
Schmidt transition!^ 

The transition rate discussed above is a dynamical 
property of the model. However, the Fermi edge singu- 
larity also manifests in ground state properties. Denote 
the ground state of the full Hamiltonian by |GS). We 
define the probability to find the impurity in state 

a) when the system is in |GS) 

n„ = (GS|(|a)(a|)|GS). (2.14) 

Let A be the index such that \F\) ® |A) is the ground 
state of i?o, the Hamiltonian in the absence of impurity 
tunnelling. It is straight-forward to show that, to second 
order in iJ^, Ua with a ^ A, is given by 

r° du: Wax{i^) 



! ( ^2- (2-15) 

27r (uiax - L^Y 

(See Appendix [X] for a derivation.) Note that, since A 
refers to the ground state of Hq, ujax is positive and the 
denominator of the integrand is therefore always non- 
zero. Using the expression of Eq. l2.11l for Wa\ one finds 



riaiuJax) = 1 7a A I 
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Here D^x is a constant of order one (at least for A^_^ < 
2). The regime of validity of the above equation is the 
same as that of Eq. l2.11l for WaX- In particular, for 
A^_^ < 2, it becomes invalid at sufficiently small ujax 
due to a break-down in perturbation theory in Ht- For 
A^_^ > 2 there are additional contributions that are an- 
alytic at small uJaX- These contributions are however 
not important for the quantities that we consider in this 
work, as they don't give rise to any singular dependence 
on LOax- For large tOaX there is a cross-over from the power 
law of Eg. 12.161 to an inverse square law Ua ~ (woa)"^- 

III. CURRENT AND FRIEDEL SUM RULE 

Using the commutation relations between Fourier com- 
ponents of the density operator, the continuity equation 
i[p{x),H] = dxj{x) is straight-forwardly derived where 
the current operator is 

j{^x)=p{x)+Y,^-^\a){a\. (3.1) 

a 

As a preliminary to our goal of studying correlation 
functions, we consider the average current when the sys- 
tem is in an arbitrary stationary state. Let D be any 
density matrix that describes a stationary state of the 
system, i.e. [H,D] ~ 0. The expectation value of the 
commutator of any operator A with the Hamiltonian is 
zero, i.e. 

{[A,H])=iy{D[A,H]} = 0. (3.2) 

From the continuity equation it then follows that 
dx = 0. Integrating from — oo to a; and using the 

explicit expression of Eg. 13. II for j{x), we obtain 

(p(x)}=p-5]^n„, (3.3) 

where Uq = {\a) {a\) is the probability to find the impu- 
rity in the state \a) and p is the density at infinity, which 
is equal to the homogeneous density in the absence of the 
impurity. The total charge AiV displaced by the impu- 
rity, is obtained by subtracting p and integrating over x. 
In terms of the scattering phase shifts defined in Eg. l2.3[ 
we obtain 

1 

AN^-Y^n^S^. (3.4) 

a=l 

For the case of a static impurity, i.e. M = 1 and 
ni = 1, the result reduces to Friedel's well-known sum 
rule AA'^ = S/tt^ Eg. 13. 41 represents a generalization of 
Friedel's sum rule to our dynamic impurity model. That 
such a generalization exists is not unexpected. A classic 
result in many-body theory is the proof by Langreth^ 
that AN = S/tt holds for another dynamic impurity sys- 
tem, namely the Anderson impurity model. 



Before considering correlation functions, we define the 
current-fluctuation operator Sj{x) through its Fourier 

transform 6jq — J^j^^^dx e^^''^6j{x), where 

a 

jaq^ Pq + ^- (3.5) 

This definition is such that the expectation value of Sjq 
with respect to the non- interacting ground state |Fx)® |A) 
is zero. From the continuity equation follows that the 
expectation value of Sjq with respect to the interacting 
ground state |GS) is also zero, unless q — 0. Using the 
definition of Aax in Eq. l2.11l and the fact that J^a ~ 
1, one furthermore finds 

{GS\Sjq\GS)^5oqJ2^o.xnc,. (3.6) 

a 

For a finite interaction strength between the impurity 
and the Fermi gas, the relative phase shifts A^a are fi- 
nite, and hence in the thermodynamic limit L — > oo, 
the Fourier transformed fluctuations operator 6j{x) has 
a zero expectation with respect to the full ground state 

(GS| Sjix) |GS) = jYl ^ Oli-oo ■ (3-7) 

a 



IV. CURRENT-CURRENT CORRELATORS 

The similarity of the generalized Friedel sum rule of 
Eq l3.4l to the static impurity result, A A = S/n, begs 
the following question: Given the probabilities Ua, can 
the effect of the impurity on the electron gas perhaps 
be accounted for by a static effective potential VcS = 
J2a''^o'Va{x) in stead of a dynamic impurity? In this 
section, we show that the answer is "No" . We show this 
by examening the ground state current-current correlator 
(GS| Sj{x)Sj{y) |GS), which can in principle be extracted 
from the statistics of repeated experiments in which the 
current is measured. (See Appendix IbI) 



A. Zero'th order 

We write 

{GS\Sj{x)Sj{y)\GS) ^ Co{x,y) + C2{x,y) + O{j*), (4.1) 

where c„ is of order of rt in 7. The zero-order term 

coix,y) = {Fx\5j{x)Sjiy)\Fx), (4.2) 

becomes singular when x — ^ y. To deal with this sin- 
gularity, it is convenient to deflne a regularised current 
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operator 

6jr{x) 



q 



_L/2 71" {x-x')'^+r 



^KT-Z2^J(^)' (4-3) 



the last line being valid for |x|/L ^ 1. At the end of 
the calculation the limit r — > 0+ is taken and results 
independent of the regularization is obtained as long as 
X ^ y. From Eas. l2.5l and 12.91 follows that 



{Fx\6jp6jg\Fx)^9ip)Sp..g 



Lp 
2tt' 



(4.4) 



Performing the inverse Fourier transform we then find 
{Fx\5j{x)8j{y)\Fx) 



—I 
2^ 



2itL 



p>0 



9. [l 



(4.5) 



Taking the limits L oo and r — 0+ we find 



co{x,y) 



1 1 
(27r)2 {y - x) 



(4.6) 



This result is independent of vx{x) and represents the 
exact result for a static impurity. 



B. Secdond order 

We obtain the second order correction C2{x,y) to the 
current-current correlator by expanding |GS) in Eq. l4.1l to 
second order in 7. The perturbation expansion of eigen- 
states can be obtained from a perturbation expansion of 
the interaction picture time evolution operator describ- 
ing the situation in which the perturbation is switched on 
adiabatically [Gell-mann Low theorem]. Thus, to second 
order in the 7's 



|GS) 

IV^ 0:10:2 ) 



\Fx) (g)\X)+J2 \^^) ® I") + XI l'/'«i"2) «) |a2) 

a oi 0:2 

-z^T^A / dte'^'Q^xit)\Fx), 



(4.7) 
(4.8) 
(4.9) 



Here 77 is a small positive constant and the limit 77 — > 0^ 
must be taken after the expectation value in Eg. 14. II is 
evaluated. The operators Qap{t) are defined as 

Q^p{t) = e'"-'e-'"f\ (4.10) 

We also define the expectation value 

P^p{t) = {Fx\Qo.p{t)\Fx) 

= e^("°^-"'3^)* {Fx\ e'^'-'e-'^^* \Fx) . (4.11) 

It is convenient to perform a Fourier transform 



L/2 

C2{p,q) = / dx 

-L/2 



L/2 



L/2 



dye-'P''e-"iyc2{x,y). (4.12) 



Expanding the ground state |GS) in 7 as in Eg. 14. 71 we 
obtain 

C2(p,g) =Xc2(p,(7), 

O 

+ {MSjxp5jx<i\Fx}. (4.13) 



We now substitute \ipa) and \i^ap) from Eqs. l4.8l and 
l4Jl into Eq. l4T3l The individual terms in Eq.|4j3 
are evaluated by exploiting the bosonic commuators 
[Sjaq,Sjpq>] = Lq5q-q>/2TT obcycd by the jaq operators, 
(cf. Eq. i2.5p . For g > 0, the 5jaq correspond to bosonic 
annihilation operators, with corresponding creation oper- 
ators Sjl^q = Sja-q- The State \Fx) is a coherent state, i.e. 
an eigenstate of the annihilation operators, (cf. Eq.[! 



Sj^q \Fx) = A^^ \Fx) , 



ZTT 



(4.14) 



The commutator of Sjaq with Qp^it) is easily calculated 
using the BCH identity. The result is 



[6ja.q.Qpj] = {e''' -l)A 



(4.15) 



After some algebra, we obtain the expression C2{p,q) = 
C^{p,q) + Ci{^p,^qr where 
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C^{p,q)= lim |7.ApALAL / dh 

jj->0+ 



X e 



»)(tl+t2) 



We change integration variables to t = ti and r = ^2 — and perform the t integral to obtain 



(4.16) 



C^ip,q)^ lim 



1 7a A I 



-A'' A'' 

riTo+ 2r] + i{p + q) 



dT [1 + {e'P- - 1) + 0(g) (e'«- - 1)] e"-P^xir). 



From Eg. 14. Ill follows that Pa\ (t) is related to the transition rate Wa\ by a Fourier transform. 



(4.17) 



(4.18) 



Substituting for \ja\\'^ PaxiT) from Eq. l4.18| performing the r integral and taking the 77 — > 0+ limit, we find 



and hence 



2^ 



T^aA(^) 



dip) 



LO - LOaX -p UJ - UJaX 



0{q) 



1 



1 



dw 
2^ 



W^aA(w) <^ sign(p) 



W - WaA - bl W - WaA 



+ sign(g) 



W - UaX - q ^ - l^aX 

(4.19) 



1 



UJ-U;aX~\q\ LO-UJaX 



Fourier transforming back to C2{x,y) we obtain 
C2(a;,J/) = 



(4.20) 



^0 

^Wax{i^) dz I dz'[Ac,x{x + z)Aax{y + z') + Aax{x + z')Aax{y + z)]ha{z'-z), (4.21) 



27r 



r 



where 

AaA(a;) 
and 

ha{z) = 

For large 



Zn 



— sin{qz) 



Vajx) - Vxjx) 

2tt 



(4.22) 



1 



^41 



LOaX - LOaX - UJ + \q\ 

(4.23) 



ha{z) ^ - 



TT {ujax - a;)3z3 



1 



z expansion of ha{z) (cf. Eq. l4.24p . For \z\ :s> x+ — x_, 
ha (z) is a slowly varying function on the scale of a;+ — x_ . 
The leading order behavior in \x — y\ of C2 (x,y) can be 
obtained by evaluating ha{z' — z) at \x — y\ 

c^(.,,).-i(A„A)^^a... r ^-i^^, 

(4.25) 

where A^a = ^2; AQ,A(a;). Substituting into the above 
the result of Eg. 12. 151 where the rate WaX was related 
to the occupation probability Ua of impurity level a we 
obtain the simple final result 



TT (ulaX 



•=2(2;, 2/) 



-(A, 



- Loy^z'^ 

(4.24) 

Since Aq,a(x) is zero for x ^ I, i.e. outside the scat- 
tering region, the first term in the square brackets in '^'-'^ min{a;, y} ^ a;_ and max{x, j/} ^ x+ 



\x - y\ 



rdu^^riaiujax), (4.26) 



Eg. 14.211 gives a non-zero contribution only when x > x^ 
and y < x+. The second term in the square brackets on 



For ujax not too large, i.e. not too far from the Fermi 
edge threshold, the energy dependence of the correlator 



the other hand only gives a non-zero contribution when can be obtained from the expression (Eq.IHH) for n„ 
X < x+ and y > x^. Thus for x and y both to the left that we discussed in the section on the Fermi edge sin- 



(< x_) or both to the right (> x^) of the scattering re- 
gion, €2(2;, y) is zero. It can be shown that this statement 
is true to all orders in Ht- 



gularity. This leads to 



(4.27) 



For X and y on opposite sides of the scattering region. For A^^-^ < 3 the result diverges close to the Fermi edge 
and such that |x — ?/| 3> x+ — a;_, we can use the large threshold. As with the rate Wax and the occupation 
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probability tIq, the divergence signals a breakdown of 
our expansion in 7. Understanding this regime is an open 
problem that we are currently working on. 

V. PAIR-CORRELATION FUNCTION 

Although it can be measured directly (see Appendix 
IB|1 . the current-current correlator that we studied in the 
previous subsection is a rather abstract quantity. In this 
subsection we therefore relate it to another correlator, the 
pair correlation function, that has an appealingly straight 
forward interpretation. The pair correlation function is 
defined as 

gix,y) = {GS\^P^iy)^PHx)i:{x)i:{y)\GS) . (5.1) 

It measures the likelihood of finding an electron at x to- 
gether with another electron at y. For x ^ y it can also 
be written as g{x, y) = (GS| p{x)p{y) |GS). Since we have 
j{x) = p{x) outside the scattering region (see Eq. l3.1|) . 
our result for the current-current correlator can be used 
to obtain a formula for g[x, y) away from the scatterer, 
namely 

(5.2) 

The negative sign in front of the second term is due to 
Fermi statistics: Given that there is an electron at x, the 
likelihood of finding another electron at nearby position 
y is less than for uncorrelated distinguishable particles. 
The term C2(x, y), as well as the higher order corrections, 
are zero if x and y refer to points on the same side of the 
scatterer, i.e. if x, y < x_ or x, y > x+. 

For max{a;; y} ^ x^ and min{a;; y} ^ X-, the approx- 
imate result of Eg. 14. 261 may be used for C2 . The sign of 
C2{x,y) is determined by —duina{Lu)\^^^ ^. It is intu- 
itively plausible that the probability Ua is a decreasing 
function oiiOax: the higher the excitation energy, the less 
likely it is to find the state |a) occupied. Refering back to 
Eq. 12. 161 in our discussion of the Fermi edge singularity, 
we see that this is indeed the case not too far from the 
Fermi edge threshold, provided A^^ < 2. Exact results 
for specific model interactions indicate that remains 
a decreasing function of ujax also when uja\ is further 
away from the threshold or Aa\ > 2. This leads to the 
conclusion that d^ix^y) is always positive. Thus in the 
presence of the impurity, the pair correlation function is 
larger than in the absence of the impurity. This a truly 
two-particle correlation effect. It cannot be accounted for 
by an increase in density around the impurity, for the fol- 
lowing reasons. The increase in pair correlations occurs 
outside the scattering region, where (GS| p{x) |GS} is un- 
affected by the presence of the impurity. Also, the pair 
correlation function increases regardless of whether the 
impurity attracts or repel electrons. Finally, the increase 
in the pair correlation function depends on distance x — y 
rather than on the distance from x or y to the impurity. 



VI. CONCLUSIONS 

We studied a model in which a dynamical quantum 
impurity is coupled to a degenerate electron gas. The 
probability Ua find the impurity in excited state \a) if 
the system as a whole is in the ground state, has an 
energy dependence (cf. Eq. l2.16p 

no. ^ (WaA)"^'^"', (6.1) 

where Wqa is the internal energy of impurity state \a) 
minus the threshold energy of the Fermi edge singularity, 
and Aqa is the relative scattering phase shift that mea- 
sures the strength of the interaction between impurity 
state I a) and the Fermi gas. The power- law in Eq. l6.1l 
is due to the Fermi edge singularity, a phenomenon that 
is usually discussed in the context of impurity transition 
rates. We note that the divergence of Ua at the Fermi 
edge threshold for A^a < 2 signals a breakdown in per- 
turbation theory in the tunnelling matrix elements 
of the impurity. 

Eq. l6.1l establishes that the Fermi edge singularity 
manifests itself in ground state properties of the impurity. 
How does it manifest itself in properties of the electron 
gas? We found that, for any stationary state, the average 
density of the electron gas is (cf. Eq. l3.3p 

M^))-/^-E^"- (6-2) 

a 

where p is the homogeneous density in the absence of 
the impurity and Va {x) is the electrostatic potential that 
the electron gas is subjected to when the impurity is in 
state I a). This leads to a generalized Friedel sum rule 
(cf. Eq. l3.4|) AA^ = naSa/ir for the number of parti- 
cles displaced by the impurity, where 5a is the scattering 
phase shift associated with Va{x). Two important points 
about these results are that (1) the effect of the Fermi 
edge singularity is contained in the occupation probabil- 
ities Ua] and (2) the average density is unaffected by the 
impurity in regions where the potentails Vaix) are zero, 
i.e. outside the scattering region. 

Having investigated the expectation values of single 
particle observables, the next logical step is to look at 
correlation functions. We therefore calculated the static 
current-current correlator (GS| Sj{x)6j{y) |GS} to second 
order in the impurity tunneling amplitudes "fap- The 
zero-order result, which is also the full answer in the case 
of a static impurity) is (cf. Ea. l4.6[) 

We obtained an expression for the second order correction 
C2(x, y) to this result fEa. l4.2l| ) that is valid for all x and 
y. It is zero for x and y both to the left or both to 
the right of the scattering region. Taking x and y on 
different sides of the scattering region, and the distance 
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\x — y\ large, we obtained the compact formula 



C2{X,V) 



Tr\x — y 



^na{uJa\). (6.4) 



Thus, current-current correlations induced by the impu- 
rity show a slow (power-law) decay as a function of dis- 
tance. These correlations are sub-leading: at large dis- 
tances, correlations that are present also in the absence 
of the impurity decay more slowly (second vs. third in- 
verse power of distance). However, the impurity-induced 
correlations can be detected by varying the impurity pa- 
rameters, as this leaves the leading order correlations un- 
affected. Due to the appearance of tt-q in Eg. 16. 41 the 
impurity induced correlations have a power law energy 
dependence ~ (wqa)^"^"^- Thus correlations are also 
sensitive to the Fermi edge singularity. The divergence 
at the Fermi edge threshold for A^^^ < 3 again signals the 
breakdown of perturbation theory in 70,^. Understanding 
correlations in the regime (Aq.^)^ < 3 and uja\ — )> is an 
open problem we are currently working on. In the regime 
where the Fermi edge singularity result for Ua (Eq. l6.1|) is 
valid, the correction C2{x,y) is positive. We have argued 
that Ha is always a decreasing function of LUaX- Thus we 
always expect a positive correction C2{x,y). 

Outside the scattering region, the current-current cor- 
relator is simply related to the pair correlation function 
gix,y)^(GS\pix)piy)\GS} as 



g{x,y)=p^-{GS\6j{x)Sj{y)\GS) 



(6.5) 



The pair correlation function measures the likelihood to 
find an electron at x together with another electron at 
y. The positivity of C2{x,y) implies that given an elec- 
tron at X, the likelihood of finding another one at y in 
the presence of the impurity increases by an amount pro- 
portional to l/\x — y\^, for X and y on opposite sides of 
the scattering region, compared to when the impurity is 
absent. It is important to note that the increase occurs 
outside the scattering region, where the average density 
of electrons is unaffected by the impurity. 



Appendix A: The relation between impurity 
transition rates and ground state occupation 
probabilities 



By changing integration variables to T = {t + 1')/2 and 
T — t — t' , performing the T-integral we then obtain 



/oo 
dre-''l^lp„A(T). (A2) 
-00 



Ua = lim 

r)-s-0+ 2r] 

Substituting for \^ax\'^ Pa\{t) from (Eq. l4.18p allows us 
to perform the r-integral. Taking the 77 — > 0^ limit, we 
obtain 



du: 

-00 27r {uJax - w)2 



(A3) 



Appendix B: Measurability of the current-current 
correlator 

We present here what we find to be conceptually the 
simplest scheme to measure the current-current correla- 
tor. Other schemes may be more general or more prac- 
tical. The Hamiltonian of Ea. l2.2l describes right- moving 
electrons propagating along the real line. A straight- 
forward realization of such a conductor is a quantum 
Hall edge state. However, another realization is a semi- 
infinite wire running from a; = — (X)toa; = 0, containing 
both left- and right moving electrons. In this realiza- 
tion, the fermion creation operators ip''{x) and ip''{—x) 
both create an electron at a position — while the sign 
of x determines whether a left-mover (positive sign) or 
a right-mover (negative sign) is created. The physically 
measurable current at — |a;| corresponds to the operator 



j(-N) = j(-N)-j(N), 



(Bl) 



where j(x) is the current operator defined in Eg. 13. II In 
other words, the physically measured current equals the 
current of right-movers minus the current of left-movers. 
In spectral representation, we write J(— |a;|) as 



Ji~\x\) = X! "^"^ 1"*^ ^"^1 



(B2) 



where J„i and \m) are the many-body eigenvalues and 
eigenstates of J(— |a;|). 

By performing current measurements at — \x\ on a suffi- 
cient number of identical systems prepared in the ground 
state, the probability density 



P{J)=Y,SiJrn-J)\{GS\ 



(B3) 



In this appendix we derive a relation (Eg. 12. 12]) be- 
tween the ground state occupation probability ria and 
the transition rate Wa\ of the impurity. The starting 
point is the perturbation expansion of Egs. l4.7l and 14.81 
for the ground state. From these equations follow that, 
to second order in the 7's, tIq,, with a ^ A, is given by 



fim |7„Ap /° dt f dt' e^('+''^P^x{t~t'). 

-00 J —00 

(Al) 



for an outcome J can approximately be determined. 
Given P{J), its second moment 



/oo 
dJ J^P{J), 
-00 



(B4) 



can be extracted. From the definitions of J(— jxl) and 
P{J) follows that 

M,Hx\)^{ji\x\f) + {jHx\r) 

-OX-N)j(N))-(j(N)j(-NI))- (B5) 



9 



In the above expression, expectation values are with re- 
spect to the ground state. 

The expectation value turns out to be position- 

independent for X outside the range of the impurity in- 
teraction, i.e. X ^l. The proof is as follows: 

d,{]{xf) = {d.Kx)j{x)) + {j{x)d,j{x)) 

= i{[p{x),H]j{x))+i{j{x)[p{x),H]) 

= i{[3{x),H]j{x))+i{3{x)[3{x),H]) 

= 0. (B6) 

In the second line we exploited the continuity equation, 
while in the third line we used the fact that p{x) — j{x) 
(cf. Eq. 13. ip for a; ^l. The last line is obtained by noting 
that the expectation value is with respect to |GS) which 
is an eigenstate of H . 



Thus, M2(— Ixj) can be measured, and is given by 

Ah{-\x\)=M+{j{~\xM\x\)) + {3{\xM-\x\)), (B7) 

where M is independent of x. Our main result (Eq. 14.26)) 
then translates into a prediction for the x dependent part 
of Ma: 

M2{-\x\)-M = -_^^(A,A)'9<.„,n„(w„A). (B8) 
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